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I. INTRODUCTION 



The Material Balance Method of calculating the behavior of pet- 
roleum reservoirs is one of the petroleum engineer's most important 
tools. Introduced in 1936, the Material Balance has undergone many 
improvements intended to simplify the mathematical techniques of handling 
the complex physical relationships involved. In spite of all advances 
in technique, the method remains unwieldy and laborious „ 

The basic concept of the Material Balance applies the Law of 
Conservation of Mass to the process of producing petroleum from an under- 
ground reservoir. The amount of hydrocarbons originally present in the 
reservoir must equal the amount produced plus the amount remaining in the 
ground at any given time during the production life of the reservoir. As 
normally used, the Material Balance is -written volume trie ally. It relates 
the volumes of oil and gas (produced and remaining) by means of their 
physical characteristics (gas volume, gas solubility and oil shrinkage) 
at the pressure prevailing in the reservoir. 

The Material Balance has been most widely used for the solution of 
two general problems t the estimation of original oil in place and the 
prediction of future production performance of a reservoir. The reservoir 
discussed in this paper will be of the Solution Gas Drive Type, with fixed 
reservoir volume, and having no gas cap and no water encroachment into 
the reservoir. 



1 



2 



A. Estimation of the Original Oil in Place 

The quantity of oil and gas originally present in a reservoir may 
be estimated by observing the drop in reservoir pressure during the pro- 
duction of a certain amount of fluids, and relating these data to the 
laboratory-measured behavior of the reservoir fluids caused by pressure 
changes. The volume of reservoir fluids (with known shrinkage and 
solubility characteristics) which will undergo a certain pressure change 
because of the removal of a given quantity of these fluids is a fixed 
quantity. Thus, the drop in reservoir pressure caused by a given amount 
of oil and gas being produced defines the volume of these fluids origi- 
nally present. 
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Bo Prediction of the Future Production 
Performance of a Reservoir 

For any given pressure assumed to prevail in the reservoir at some 
future time, the volumes of oil and gas produced (v/hen related to their 
physical characteristics at that pressure) constitute a unique solution 
to the equation balancing the original volume of fluid with the volumes 
produced and remaining. An additional concept is involved however. 

The permeability of the porous reservoir medium to oil and to gas changes 
as production proceeds. Thus, as the saturation of liquid in the reser- 
voir decreases as oil is produced, the relative amounts of oil and gas 
which flow to the well will change. This phenomenon cannot be directly 
related to the Material Balance itself, yet it determines the relative 
quantities of oil and gas which will be produced. The relationship be- 
tween the relative permeability to oil and gas and the liquid saturation 
in the reservoir must be known or assumed, since it must be considered in 
arriving at the amounts of oil and gas to be produced. These amounts of 
oil and gas can then be used to predict a solution to the Material Balance 
at some assumed future point of pressure and production. 

It will be the purpose of this thesis to develop a graphical means 
for predicting the performance of a depletion drive reservoir based on a 
modification of the Material Balance equation. 



II. REVIEW OF THE LITERATURE 



The use of the Material Balance for estimating oil in place was 
developed by Schilthuis in 1936. 1 His equation (shown here for a solution 
gas drive reservoir without gas cap, water encroachment or water production) 

* 



N 



n[u + (r n - s Q )v ] 
u - u Q 



( 1 ) 



still remains the most used basic form of the Material Balance. 

O 

At about the same time Katz proposed a tabular method for evaluat- 
ing oil in place. This method was later shown by Pirson^ to be equivalent 
to the Schilthuis method. 

The use of the Material Balance Equation for prediction of future 
reservoir performance was proposed considerably later. Babson^ developed 
a trial and error solution based on the Material Balance Equation and the 
Instantaneous Gas-Oil-Ratio Equation. His method was very cumbersome and 
required a great deal of computation. It is rarely used today. In the same 
year, 19hh, Tarner proposed a solution to the problem which was consider- 
ably simpler to use. 

Muskat^ developed a form of the Material Balance expressed as a dif- 
ferential equation and applied it to prediction of depletion drive reser- 
voir performance. 

The most widely used of the Material Balance prediction techniques 
is the method using a trial and error solution of Schilthui^ Equation^* 

This method requires the simultaneous satisfaction of three equations ; 



^"References are listed in the Bibliography. 
Nomenclature shown in Appendix I. 
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1. The Material Balance Equation 



N - 



n[u + (r n - s q )v] 



( 1 ) 



2. The Instantaneous Gas-Oil-Ratio Equation 



•"o 0 K g 



+ s 



( 2 ) 



L 




3. The Liquid Saturation Equation 



S L = % + t 1 " S ff) ^ 




(3) 



manners 

1 . 

2 . 

3. 

In 

5 . 



trial and error solution may be conducted in the following 

Assume a certain drop in reservoir pressure 

Estimate the volume of oil which would be produced during 

this pressure drop 

Using this estimated value of oil productions solve the 
Liquid Saturation Equation* (3) s and obtain the corresponding 

value of _jL from the known or assumed vs. Sr relation- 

K_ k L 

o A o 

ship. 

Solve Equation (2) for the value of Instantaneous Gas-Oil - 
Ratio which would obtain at the assumed pressure. 

Using the assumed and estimated values of oil production* gas- 
oil-ratio and fluid characteristic solve the Material Balance 
Equation. If a balance is obtained* the assumed values were 
correct. If the equation does not balance* a new value of oil 
production or pressure is estimated and another trial is made. 
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This process is repeated until the equations are satisfied at sufficient 
pressure points to define the future performance of the reservoir*, 

The use of the Schilthuis Equation, while more flexible in the 
handling of such conditions as water encroachment and gas cap expansion, 
is extremely tedious. For the simplest reservoirs, twenty-five columns 
of calculations must be used, and five or six trials at each of seven or 
eight pressure points is not unusual. All this must be performed with a 
calculating machine. Despite this, it is still the most widely used 
mathematical prediction technique. 

In January, 1955» a prediction method based on the Schilthuis 
method, but greatly simplified, was published by Tracy. ^ The Tracy method 
is the single biggest advance in the use of the Material Balance Equation. 
This method will be described for the case of a depletion drive reservoir 
without gas cap or water encroachment. 

Tracy took the Schilthuis Equation 

n[u + (r n - s )v] 

N - — 

u - u 0 




and substituted values for u, u^ and r as follows s 

o n 

u = ft + (s Q - s) v 
,J o ' £ o 




This substitution resulted ins 




00 

(5) 

(6) 



N - 



(?) 
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which nay be expanded to; 




Thus, n and G are multiplied by coefficients which are functions of pressure 
only. These coefficients are ; 



4 * 




(9) 




( 10 ) 



Equation (8) then reduces to: 

N - n i + G i (11) 

b 

In this form, the equation may be used to estimate original oil in place. 

Tracy's prediction method uses this equation in two forms with N 
taken as equal to one barrel of stock and oil. 

1. As the Prediction Equation; 



An 



1 " ( n i-l + G i-1 ^g) 



K + 



i +r i-l 



( 12 ) 



2 a As the Material Balance Equations 




+ G . 



1 



(13) 
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Equations (12) and (1.3) are used in conjunction with Gas-Oil~Ratio 
Equation; 



Mo 

r i = 4g ^ K” 



s 



( 2 ) 



and Liquid Saturation Equation; 

s l - s# + c 1 - v f^i^y «> 

K„ 

Equations (2) and (3) being linked by knowledge of the rr- vs. 

K o 

ship. 

The use of Tracy’s Prediction Method calls for the following 
procedure s 

1. Starting with a point in production history (with known values 
of n^* </>, r^ for a given pressure)* assume a new lower pressure point. 

2. Estimate an instantaneous gas-oil ratio for this pressure 

3. Using this pressure and gas-oil ratio in the Prediction Equation 
(12) , compute a value for An. 

U. Add this value of A n to - 1 (n^ from the previous pressure) 
to get the value of n^ for the new assumed pressure. 

5. Solve the Liquid Saturation Equation* (3)* using this value of 

n-; and obtain a value for Jg . 

K K ° 

6. Using the value* solve the Gas-Oil Ratio Equation. This 

K 

o 

value of Instantaneous Gas-Oil Ratio should equal that estimated in step 2. 
If it does not* use the value of found in step 6 as the estimated value 

in step 2 and re-do steps 2 through 6. When the values of Instantaneous 
Gas-Oil Ratio in steps 2 and 6 are equals substitute the values of n . and 
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G^ into the Material Balance Equation (13) and solve. The sum of 

n-; + G., </> should be 1.00 + 0.02 e 

i n i g — 

The number of trials required is cut down to two or three in Tracy's 
method because Instantaneous Gas-Oil Ratio - a relatively insensitive 
factor - is estimated. This cannot be done in the Schilthuis Equation 
as the cumulative value of gas-oil ratio, rather than the instantaneous 
value, is used in the equation. The cumulative gas-oil ratio value cannot 
be estimated with any degree of accuracy,, 

Sturdivant, a ser; j_ es 0 f three articles, expanded upon the 
use of the Tracy method developing a formula to be used in the event of gas 

Tj 

reinjection into the reservoir. It was also shown that accuracy to three 
significant figures in the calculation was sufficient. 



Ill, STATEMENT OF THE PROBLEM 



This investigation was conducted with the object of developing a 
graphical method for solution of the Performance Prediction problem. It was 
felt that the Prediction Equation developed by Tracy* with its separation of 
the pressure factors into coefficients* would lend itself to such treatment. 

In the Schilthuis Equation* the complexity of .form appeared to preclude any 
simple graphical treatment. 

Of the various graphical problem solving methods available* nomography 
appeared as the best for trial and error solutions. The forms of the equa- 
tions indicated that nomographs of reasonable simplicity could be constructed 
to solve them. For trial and error solutions* it was felt that the norao- 
graphs should be simple enough that the operator could see at a glance how 
much a given change in one factor would affect the other factors in the prob- 
lem. Also* with simple* one-setting nomographs* the operator can readily 
enter the diagram with the dependent variable (or answer) and work backwards 
to find the value of independent variable (or normal input). This is of value 
on trial and error solutions. For these reasons* as well as for ease of con- 
struction* it was desirable that the nomographs be of the type in which a 
single line crosses all scales* obtaining the answer in a single setting. 

There are a number of nomographs in the literature concerned with 

~1'3 1 c* ~i z' 

various aspects of the Material Balance problem 5 . All of them are 

quite complex* requiring entry with several individual pressure variables* 
and requiring several settings of a line before a solution is obtained. 

These nomographs are general in application and could be used on reservoirs 
having the same ranges of variables as the nomographs. 
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By constructing nomographs for a specific reservoir* all of the fluid 
characteristics (which are functions of reservoir pressure) and combinations 
of them* can be treated as a singles or recurring* variable of pressure. 

This reduces the number of variables to the point where simple., single^ 
setting nomographs can be constructed. Construction of the nomographs for 
a single reservoir allows selection of optimum scale lengths to specially 
handle the ranges of the variables for that reservoir. More generalized 
nomographs would require the ranges of the scales to be broad enough to fit 
all or several reservoirs. This would result in less accuracy for any 
particular reservoir. 



IV. METHOD OF INVESTIGATION 



The four equations used in the Tracy method of prediction were 
separately investigated to find the forms most likely to yield the simplest 
nomographs. These equations;, given previously as Equations (12) , ( 13 ), (2) 
and (3) -will be considered separately. 
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A. The Prediction Equation 



17 

This equation was given as g 



An = 



1 * < n i-l K + Vl O 



where 



4 * r a 



r-i + r 



i T *i-l 



(Hi) 



(15) 



It is readily rearranged into this form; 



(n. + A n) f! + (G. , + r A n) p = 1 

v i“l n v i-l a p 



( 16 ) 



In this form a it is seen to be equivalent to the Material Balance 
Equation as s 



and 



+ An s n. 



G i-1 r a A n = G i 



(17) 



( 18 ) 



If these values are substituted into Equation (16) 9 the result is the 
Material Balance Equations 

n. ^ + G„ 0 = 1 ( 13 ) 

in i g 

Returning to the rearranged form of the Prediction Equations 

(n. , A n) fS + (G + r A. n) 6 = 1 (16) 

l-l n l-l a g v 

it is noted that there are four variables s 



Ik 



(ru ^ + An), a function of oil production 
Let this be f., (u) 

(Gi_i + r a An), a function of the gas produced with the above 

amount of oil 
Let this be (v) 

^ , a function of reservoir pressure 

Let this be f^ (W) 

</> , a different function of reservoir pressure 

s 

Let this be (W) 

Thus, Equation (16) can be written as s 



i\ (u)x f 3 (W) + f 2 (v) * f u (W) - 1 



(17) 



Note that the pressure, represented here by tT W" is a recurrent 
variable . 

Divide Equation (17) by f (W) , thens 



% (W) q 

f ! (u) + f 2 (v) ' fTW ■ f~TWJ 



( 18 ) 



Nov/ let 



f 3 (w) 



= f^ (W), a different function of pressure 



and 



q <vt 



= f^ (W), another function of pressure. 



Equation (18) now becomes 



(v)*f£ (W) = f 6 (W) 



q (u) + q 



(19) 
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This form is capable, of representation by a simple nomograph having 
the following configuration? 



n scale 

( n i “ n j_i + An ) 



G scale 

(G. = G. + A n r ) 
x i i-1 a' 




The geometric proof for this nomograph is given in Appendix IV. 

Now, evaluating the pressure functions f ^ (W) and f (W) results in 

? 6 



(W) 



(W) 

r 3 Tiy 



4 

Jl 

V 



( 20 ) 



l v 



sA 



f/c w * r - i 






\ v 



0 

S I 

y 



(21) 



These values, when substituted into the Prediction Equation, gives 



( n i~l + & n ) * 



< G i-i + r a ^ n) 



C4 



) 



s) 



f (&0 _ 

1 - l V -o 

(4 - •) 



( 22 ) 



Equation (22) is the form of Tracy® s Prediction Equation which is 



solved by the Prediction Equation Nomograph 
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Bo The Material Balance Equation 



This equation is equivalent to the Prediction Equation and can use 
the same nomograph 0 One solution of the nomograph serves to satisfy both 
the Prediction Equation and the Material Balance Equation „ 

Placed in the form used to construct, the nomograph, the Material 
Balance Equation becomes s 



n „ + 

1 




Equation (23) is solved by the Prediction Equation nomograph,, 
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C e The Liquid Saturation Equation 



The Liquid Saturation Equation is used to obtain a value of ^ from 

K 0 

a known and extrapolated relationship between S T and a 

Ij K c 



In the equation 



S * S+ (1 - SL) 
L iff W 



) 



V^o 



( 3 ) 



The values of S and ft a remain constant » 

W 

Let Sy be K x and let be K, 



The equation can then be rearranged: 



S L - Kj - (K 2 ) (/?) (1 - np (21.) 

or 

Log - K l ) - Log K 2 * Log (1-np (25) 



Equation (2A) was used to construct a parallel line nomograph with 
logarithmic scales having the following f orm 

Log (1-n^ Log (Sj-K^ Log (K 2 P) 



i 



-F or proof of the geometry of this nomograph 5 see Appendix IV. 
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Log (l~n^) is a simple function of n^ and the left scale was therefore 

graduated directly in units of n^. Similarly* {3 , and therefore Log K 2 ^ 

is a function of pressure, so that the light scale was graduated in units of 

pressure. The center scale was graduated on one side in units of and by 

use of the known _£ vs. relationship, the other side was graduated in 

K_ 

K ® 

units of . The finished nomogram appeared as follows: 

K o 

-P 






<•3 /? 5 ' 



alnter 



i 



’''Read 



( 'Enter 
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D. The Gas-Oil Ratio Equation 



In the instantaneous Gas-Oil Ratio Equation 



A-0 

Al% V fa 



*+• 3 



( 2 ) 



/ ^° , ft , v» and s are specific functions of the reservoir pressure. 
By combining all of these but s into a single function of pressures the 
equation was rearranged to; 



A- 

i x(F) = 



( 26 ) 



This form is identical to 



f x (u) - f 2 (v) . 1* 3 (W) * f (7T) 



( 27 ) 



where f (W) and f^ (W) are different functions of the recurrent variables 
reservoir pressure. 

A nomographs similar in theory ^ to the Prediction Equation nomograph 
was constructed „ The finished form appeared as follows; 




♦See Appendix IV for geometric proof 









\ 





20 



E. Average Gas-Oil Ratio Equation 



It was found desirable to construct a very simple nomograph to solve 
the equation 



r, + 
1 



r i-l 



2 



( 15 ) 



The equation was rearranged to? 



r + r = 2 r 

i i-1 a 



A nomograph with uniform parallel scales set equal distances apart 
was constructed. 



Cc 






4 - 



/ 




Entry can be made on any two scales and a value read on the third 



scale 



V. USE OF THE NOMOGRAPHS TO PREDICT THE PERFORMANCE 



OF A SPECIFIC RESERVOIR 

To test the application of the nomographic prediction method* an 
actual reservoir was selected and nomographs constructed using the data 
available concerning this reservoir. The particular reservoir chosen was 
a mid-continent limestone reservoir* which is still in the early stages of 
production. For this reason, it will be called Reservoir r, X. n The basic 
data available on Reservoir :, X‘* are listed in Appendix V. In addition* two 
sets of prediction curves (Pressure and Gas-Oil Ratio Vs. Cumulative 
Production) were available for comparison with the results obtained in this 
investigation by the nomographic method. One of these sets of prediction 
curves was obtained by trial and error use of the Schilthuis Equation. The 
other was obtained by using Tracy's method. 

The insensitivity of instantaneous gas-oil ratio values over a range 
of incremental oil production values was mentioned in connection with the 
Tracy method. The relative positions of the n, G and P scales on the 
Prediction Nomograph show this feature graphically. The pressure scale is 
very close to the gas production scale. For this reason, it is best to enter 
the nomogram with an assumed value of oil production, rather than to assume 
a gas-oil ratio as in the tabulated Tracy method. 

The following approach was found to give the easiest solution using 
the nomograms : 

lo Set up a tabulation sheet as follows % 



21 



22 



(0) 


P 


1209 

(A known point of 
Production History) 


1100 


900 


(1) 


A n 




0.00773 


0.0121*8 


(2) 


n i-l 




0.0319 


0.03963 


(3) 


n i 


0.0319 


0.03963 


0.05211 


(M 


G i-1 + r a 


1)6.75 


65.5 


108.3 


(5) 


G i-1 




16.75 


65.5 


(6) 


r a An 




18.75 


1)2.8 


(7) 


r a 




2h25 


3l)30 


(3) 


r i-l 




21*00 


21*50 


(9) 


r i 


21*00 


21*50 


hh20 


( 10 ) 


h 

K 0 




0.0362 


0.085 


(ID 


r i 


(From G.O.R. 
Equation) 


21*50 


14^20 


(12) 


(n ± ) (1.8 x 10 7 ) 




713 o0 


937.5 


(13) 


+575,000 




575.0 


575.0 


(IK) 


Cumulative Production 
(in 103 bbls.) 




1288.0 


1512.5 



2. Select a point of actual production history, knowing pressure, 
oil production from the bubble point to that pressure, n^, expressed as a 
fraction of original oil in place, gas production from the bubble point to 
that pressure expressed as a function of original oil in place, and instan- 
taneous gas-oil ratio, r^. 

3. Assume a new pressure, 100-200 psi lower than the historical 
point used. The Hj_, and of the historical pressure point now become 

n i~l, ^i-l* r i 1 an a b° ve tabulation. 

1*. Estimate an increment of oil production, An. Enter this on 
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line one of the tabulation. Add to it the n. _ . The sum is n.. 

i~l 1 

5. Using this value of n^_ at the assumed pressure, enter the pre- 
diction nomograph on the n and P scales and read the value of 

+ A n r a on ^h 0 G scale. Enter this value on line it of the 

tabular sheet and subtract the G. _ from it. The difference is An r . 

i-l a 

6. With a slide rule, divide the value of An r , (line 6) by 

cl 

An (line 1) to get r a (line 7). 

7. Enter the average gas-oil ratio nomograph with r^_^ and r a and 
find r-j_ (line 9). 

. K g 

8. Using m (line 3) and Pressure, find from the liquid 



nomograph and enter this value on line 10. 



K, 



9. Enter the gas-oil ratio nomograph v/ith _I and Pressure and 

K o 

find the r a value based on the gas-oil ratio equation. Place on line 11. 

10. Compare lines 9 and 11 which will be equal if the assumed value 
of An was correct. 

Note that as assumed values of An are increased, resulting values 
of on line 9 (from prediction nomograph) decrease sharply, whereas values 

of Tj_ on line 11 (from gas-oil ratio nomograph increase slowly. Thus, if the 
r^ of line 11 is less than that of line 9* a slight increase of A n will 
bring them together. If the two values of r^ are within a few hundred of 
each other, set the A n r a value from line 6 on the slide rule ll D u scale, 
place the value of r (line 7) opposite it on the H C l * scale, and read a 
nearly- correct value for An under the index on the f *D fl scale. 

If this new value of An causes a radical change in n^, it may be 
necessary to enter the nomograph for a new G^__^ + A n r & . However, from 
the scale positions, it can be seen that G^_ n + An r a is very insensitive 
to changes in An at most pressures. 



11. If the new value of An caused more than a slight change in 



n , re-enter the liquid saturation and gas-oil ratio nomograms to obtain a 
i 

corrected r^ for line 11* 

12. Lines 9 and 11 should now be within perhaps 5>0 or 60 of each 
other. Alter An slightly to bring line 9 in exact agreement with line 
10. Add the corrected A n to n^_ 1# The result is n^ 0 

13. Multiply n^ by the volume of original oil in place for line 

12 and add line 13 * the production from the original pressure to the bubble 
point. Line ill is the predicted cumulative production at the assumed 
pressure. 

lU. Assume a new pressure about 200 psi lower and repeat the 

process. 

Special care in the averaging of r^ and must be taken when 

computing at about the peak of the Gas-Oil Ratio Vs. Production Curve , as 
any averaging process assumes local linearity of a curve. It is best to 
use smaller pressure increments in this region. 

Using the method described above, about two hours were required 
to calculate the predicted performance using seven pressure steps. 

An alternate method was developed which takes about the same time 
to use, but which may be a little more accurate g 

1. Use the same tabulation sheet as shown on page 22 above. In 
addition, a large "scratch graph sheet" of squared paper (10 x 10 to the 
inch was used) and a sheet of scratch paper are necessary. 

2. Same as step 2 above. 

3 . Same as step 3 above. 

li. Estimate three values of An quite close together. Jot on 
scratch sheet or carry mentally. 
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5. Enter the prediction nomogram with the middle value of A n 

and read a value of G, , + An r . (It will be seen that the change in 

l-l a 

G i-1 + A n r a over the three values of A n is so small as to be un- 
readable.) 

6. Subtract G. ■, from G. _ + An r . Divide An r by each of 

i-l l-l a a J 

the 'three values of An from step h obtaining three values of r . 

cl 

7. On the scratch graph sheet, plot the three values of r vs, 

cl 

the respective An values. 

8. Add the three values of An to n^_^. Use the resulting 3 values 
of nj_ in the liquid saturation anc gas-oil ratio nomographs to obtain three 
values of rj_. 

9. Using the average gas-oil ratio nomograph and entering with 

r. n and the three r. values of step 8, find three values of r . 
i-l l r a 

10. Plot these values against A n values on the scratch graph 

sheet. 

11. Where the lines of step 7 and step 10 cross, read the correct 
values of A n and r a . Enter these values on the tabular sheet and check 
the computation on the nomographs, comparing tabulation lines nine and 11. 

12. Make minor adjustments of A n to match line 9 to line 11. 

13. Same as step 13 above, 
liu Same as step lb above. 

Figure 1 is a graph comparing the results obtained by the nomo- 
graphic method with those obtained by trial and error solution of the 
Schilthuis Equation and by Tracy’s method. The results in general are 
in good agreement. The higher gas-oil ratio values shown on the Schilthuis 
method curve are thought to have resulted from use of different values of 
”v” than were used in the Tracy method and in the construction of the 
nomographs. The Schilthuis method curve was taken from an engineering 
consultant’s report on Reservoir H X ! ' . 
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Comparison of Reservoir Performance Prediction Results 
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VI. SUUMARY AND CONCLUSIONS 



This investigation had as its object the development of a graphical 
technique for predicting production performance of a petroleum reservoir. A 
graphical technique involving the use of nanographs based on the Tracy method 
was worked out. Using this method nomographs were constructed to predict the 
performance of a specific actual reservoir called Reservoir M X." The per- 
formance of Reservoir "X" under depletion drive was predicted using these 
nomographs. The results, compared with the results obtained by other 
methods, show that the use of nomographs constructed for a specific reser- 
voir gives a forecast of reasonable accuracy. 

The total time involved in nomograph construction and in solving the 
prediction problem with these nomographs proved to be about the same as that 
required to solve a reservoir problem by the Tracy method. However, reser- 
voir performance predictions are often repeated, starting with a later point 
in production history, as the field develops. With this in mind, it is be- 
lieved that the nomographic method would save considerable calculation, A 
repeated prediction could be run in about two hours using the same nomo- 
graphs. A repeated calculation, even by the Tracy method, would take 
oenaldg*afely longer. 

As the production history of the field develops, it may be found that 
the basic field data relating relative permeability ratio and saturation were 
extrapolated in error. In this event, new and more correct values of 
may be marked on the center scale of the liquid saturation nomograph 
without affecting its aocuracy for subsequent use. 

The limited time available precluded investigating the application 
of this technique to a reservoir undergoing gas reinjection. It is almost 
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certain that a set of reasonably simple nomographs could be constructed to 
handle predictions under various conditions of gas reinjection* 

The time available was also too short to allow the method to be ap- 
plied to reservoirs of differing characteristics* In particular, it should 
be tried on a sandstone reservoir, on a reservoir with a high initial pres- 
sure, and a reservoir with a low initial pressure. 
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APPENDIX I 
Nomenclature 

sOL% v 

Cumulative gas production,, standard cubic feet 
G at reservoir pressure being considered 
G at previously considered reservoir pressure 

Relative permeability ratio, a dimensionless number 

‘i'otal original oil in place, in stock tank barrels 

Cumulative oil production, expressed either in stock tank barrels 

(or as a fraction, when N is taken as equal to one barrel) 

The value of n at the reservoir pressure being considered 
The value of n at the last previously considered reservoir 
pressure 

Average gas-oil ratio, standard cubic feet per barrel of stock 
tank oil 

Instantaneous gas-oil ratio at reservoir pressure being con- 
sidered standard cubic feet per barrel of stock tank oil 
Instantaneous gas -oil ratio at last previous reservoir pressure, 
standard cubic feet per barrel of stock tank oil 
Cumulative gas-oil ratio, or the ratio of total standard cubic 
feet of gas produced to total barrels of stock tank oil 
produced 

Total liquid saturation in the reservoir 
Oil saturation in the reservoir 



Saturation of connate water in the reservoir 
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Solution gas-oil ratio, or the solubility of gas in crude oil 
at reservoir pressure and temperature, cubic feet per barrel 
Value of s at original reservoir pressure 

Volume in the reservoir occupied by one barrel of stock tank 

oil plus all the gas originally in solution in that oil, or the 

two-phase formation volume factor 

Value of u at original reservoir pressure u q ~& q 

Barrels of free gas space occupied in the reservoir by one 

standard cubic foot of gas 

Formation volume factor, or the volume occupied at reservoir 
conditions by one barrel of stock tank oil plus its dissolved 
gas 

Value of $ at original reservoir- pressure 

Increment of oil production between two reservoir pressure 
steps 

Viscosity of reservoir oil at reservoir conditions, centipoises 
Viscosity of produced gas, at reservoir pressure and temperature 
centipoises 
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APPENDIX III 

Considerations Affecting the Accuracy of Nomographs 

1 . Dime nsions of the Diagram 

The accuracy of a nomograph can be increased by increasing the 
dimensions of the diagram, within reasonable limits. This is due to the 
increase of relative accuracy and fineness of scale graduations. 

2 . Order of Accuracy 

Nomographs, in general, have a degree of accuracy exceeding that 
of a slide rule of similar dimensions. In the nomograph, the ranges of 
the scales are tailored to the ranges of the variables involved. Short 
ranges are often expanded into long (and, therefore, more accurate) scales. 

3 . Decrease in Accuracy due to Complexity 

The accuracy of a nomograph decreases sharply when more than one 
setting of a line must be made to obtain the answer. The more inter- 
related lines that must be drawn on a given nomograph to obtain a solution, 
the less accurate the answer will be. 

U. Relative Placement of Scales 

The accuracy of a nomograph is affected by the relative positions 
of the scales. If possible, the scale on which the answer is to be read 
should be between the other two scales. 
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APPENDIX IV 

Geometric Proof for Nomographic Forms Used 

A. Prediction Equation and Material Balance Equation 

The equation form for the nomograph is: 

f (u) + f (v) y f (W. = f ('.V) (19) 

Id bo 



Take tv;o parallel axes for the variables u and v, and a 
curved axis in general for the variable 7/: 



u w 




The line AB, of length k, joins the zero points of the scales for u and v. 
Points x and y are any points on the u and v scale, x = m^ f^ (u) and 
y = m 2 fg (v), where an d ir^ are the respective scale moduli. 

Take any index line CD intersecting the scale for W at E. Con- 
struct EF parallel to the u and v axes. Let EF = z and AF = z’. From 
similar triangles ECH and EDG : 
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x-z _ z-y 
z’ k-z' 



(k-z f ) x = kz - zz' + zz’ - yz* 



(k-z') x + yz' - kz 



/ z> \/ •, / k V , 

x + - r i- v ' = • i r ‘l z : 

•.k-z'/-. J 'k-z 1 j \ J 



Substitute for x and y; 



/ k 



m i f i (i:) + (k^JK f 2 ( v >k- 'Wii 



i- (u) . p 

i V m i k-z 1 



L 



< v) z ) 



This equation -will become the original equation if f^ (W) and f^ (W) have 
the f oil owing values : 



% (W) 

f 6 W 



n^z 



m^( k-z » ) 



Solving for z and z 1 : 



z 



mg f 6 (70 

, m2 

% W 51 



Z ! 



k f^ (W) 

mo 

FOO + ~ 

p mq 
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Points on the curve for the W-scale are obtained by taking dif- 
ferent values of W and determining the corresponding values of z and z’. 
These points are plotted and a curve may be drawn through them* 
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B. Gas-Oil Ratio Equation 
The form of the equation is: 

f-L (u) - f 2 (v) . f 3 (VO = £ h 00 (27) 



This equation is similar in form to that discussed under "A Prediction 
Equation," and is subject to the same proof, except that the scale for £2 ( v ) 
runs in the opposite direction from the f^ (W) scale: 




After substitution for x and y, this reduces to: 





k 

m-j^ ( k-z • ) 



f x (W) - 



ko 



and z and z' are as before; 



*2 £ k (W) 

f, (W) + 

* “1 



i = 



k f 3 (W) 
f 3 (W) + 



H^2 

m n 



Points on the W-scale are determined by taking different values 



of W and computing the corresponding values of z and z'. 



U2 

Let ABC be a base line, perpendicular to the three scale axes, and 
let a and b be the distances between scales* 

Draw any index line in general, FH, forming similar triangles DEF 
and DGH, where DE and GH are parallel to ABC* 

From these similar triangles: 

x - z _ z - y 
a b 



which rearranges to: 



x y _ z 

a b ab/~(a+b) 



Substituting particular functions for x 3 y, and z : 



m l ^ ( m 2 ( v ) m 3 ^ W 

a b ab/(a+b) 



To reduce this to the original form, let 



m. 



b = 



ab 



a + b 



m2 

m^ , or m^ 



“l m 2 
“l + m 2 



This also defines the scale moduli in terms of each other, showing 
that after m^ and are selected, m^ may be computed from them* 

To find the position of the center scale in terms of distance between 
the outside scales, note that the distance from A to B in terms of AC is: 



a + b 



"i 

m, + m- 



(AC) 



(AC) 



or 
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D. Average Gas-Oil Ratio Equation 



The equation fonn is: 



r. + r = 2 r 

i i-1 a 



or 



(u) + f 2 (v) = f 3 (W) 



The proof is identical to that for the liquid saturation equation 
nomograph. 
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Data Concerning Reservoir n X" 

A Mid-Continent Limestone Reservoir 



Basic Data 

Average Depth 

Original Reservoir Pressure 
Reservoir Temperature 
Average Connate Water Saturation 
Original Formation Volume Factor, 0 
Original Solution Gas-Oil Ratio 



5300 ft. 



1773 psig 
128° F. 



33 »h% 



1.280 

ft 3 



500 



bbl. 



Data 








Pressure 


S, 




M 0 


(psig) 


(SCF/Bbl.STD) 


F.V.F. 


(C.P.) 


1500 


500 


1.280 


1.28 


1250 


1A8 


1.260 


1.32 


1000 


392 


1.2U0 


l.iil 


750 


336 


1.218 


1.5U 


500 


278 


1.193 


1.72 


250 


208 


1.16U 


2.0 


100 


iU5 


1.136 


2.37 


0 


0 


1.07 


3 oh 



Volume tricall y Weighted Average Bottom Hole Pressures 



U5 



BHP 



Date 


(Psig) 


Original 


1773 


10-3-52 


I651t 


U— 13— 53 


1558 


7-6-53 


1517 


10-1-53 


1473 


1-U-5U 


llilll 


7-6-514. 


13U2 


1-3-55 


1209 



Reservoir Fluid Characteristics ; 

Gravity of stock tank oil 

Saturation pressure 

Solution gas-oil ratio 

Formation volume factor 

Viscosity of reservoir oil 

Analysis of Reservoir Fluid Sample 

Methane and Lighter 

Ethane 

Propane 

Iso-butane 

N-butane 

Iso-pentane 

N-pentane 

Hexanes 

Heavier 



Cumulative Oil 
Produc tion 
(Bbls.) 

152,330 
353s 995 
bk3,999 
51*3 * 7 83 
631,761* 

866,063 

1,086,797 

1*1.5° API at 60° F. 
ll*87 PSIG at 128° F. 

U8I4 Std. Cu.Ft./Resid. Bbl. 

1.28 Res. Bbl./Resid. Bbl. 

1.28 Centipoise at ll*87 PS I + 128 °F, 



Weight % 


Mol % 


3.39 


23735^ 


1.70 


7.06 


3.75 


10.62 


0.82 


1.77 


2.55 


5.U6 


1.38 


2.38 


1.52 


2.63 


2.81 


U .07 


82.08 


39.66 



100.00 



100.00 
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Gas Viscosity 

Pressure 

(Psig) 



1300 


0.01477 


1100 


0.01401 


900 


0.01355 


700 


0.01280 


5oo 


0.01228 


300 


0.01193 


100 


0.01158 



Reservoir Gas Analysis 



Analysis of Gas Sample 


Mol % 


Methane 


68.5 


Ethane 


8.0 


Propane 


7.6 


N-Butane 


2.5 


Iso-Butane 


1.2 


N-Pentane 


0.5 


Iso-Pentane 


0.8 


Cyclo-Pentane 


0.2 


Hexanes plus 


o.4 


Nitrogen 


9.9 


Oxygen 


0.1 


Helium 


0.1 


Hydrogen 


0.1 


co 2 


0.1 


h 2 s 


0.1 

100.00 



Gas gravity at 60°F. + 14.7 psi = 0,819 (Air = 1,000) 



U7 



Relative Permeability Data 



Summary of Production: 


Year 


Producing 
Dwells at End 
of Year 


1951 


5 


1952 


26 


1953 


29 


1951 * 


31 



Production prior to bubble 



S L 


% 


{%) 


K o 


98 


0.01 


97 


0.0181 


96 


0.033 


95 


0.057 


90 


0.373 


85 


1.15 


80 


3.20 


75 


8.0 



Production 



Oil (Bbls.) Gas (MCF) 



Year 


Cumulative 


Year 


Cumula tii 


27,295 


27,295 


13,655 


13,655 


185,930 


213,225 


93,105 


106,760 


1 * 11 *, 03 I* 


627,259 


251*, 261 


361,021 


* 56,238 


1,083,1*97 


753,605 


1,111*, 626 


point : 


575*000 bbls. 


STO 





Total Volume of Original Oil in Places 

N = 18.000,000 Bbls. STO (Calculated by Material Balance and 

verified by volumetric value) 

On Basis of N = 1 Bbl. STO, 

Oil Production, B.P. to 1209 psig 
n = 0.0319 Bbl. 

Gas Prodiction, BP to 12 09 psig 
G = 1*6.75 SCF 

Instantaneous Gas-Oil Ratio at 1209 psig 
ri = 21*00 SCF/Bbl. 
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Calculations for Nomograph to Solve Gas-Oil Ratio Equation 
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Computation of Scale Lengths for Liquid Saturation Equation 



Log (S L - Kj) = 


Log (Kg 0 ) 


+ Log (1 


■~n±) 




K 1 = % = 0.33k K 2 


u-V 
' TT 


0.666 

17255 


= 0.5203 




Scale Labeled 








Range of 


Label Limits 


Variable 


Limits 


Function 


Function 


P 1300 to 




1.26k 


log Kp 


9.81803-10 


100 




1.135 


u 


9.77122-10 


12" 




0.129 




0.01*681 ' 


apart 




1-ni 


logtl-n.^ 


9.98677-10 


n 0.03 


n. 


0.97 


0.10 


1 


0.90 


9.951*2l*-10 

0.03253 




6'' 










from others S^ 


S L 




log^-I^) 






(1) 


(2) 


(3) 


(k) 




P(psig) 


0 


k 2 * 


Log 




1300 


1.26k 


0.6577 


9.81803-10 




1200 


1.256 


0.6535 


9.81525 




1100 


1.21,8 


0.61*93 


9.8121*5 




1000 


1.21*0 


0.61i5>2 


9.80969 




900 


1.232 


0.61*10 


9 0 80686 




800 


1.222 


0.6358 . 


9.80332 




700 


1.213 


0.6311 


9.80010 




600 


1.201* 


0.626U 


9.79685 




500 


1.193 


0.6207 


9.79288 




k00 


1.183 


0.6155 


9.78923 




300 


1.171 


0.6093 


9.781*83 




200 


1.156 


0.6015 


9.77921* 




100 


1.135 


0.5905 


9.77122-10 


Limiting Values of S^: 


(1) 


(2) 


(3) 


(1*) 


P = 1300; n ± = 0.03 


1-S}f 

~^7 


(1-n) 


(l) x (2) 


(3) + %=S£ 


P = 100; n ± = 0.10 


0.S203 

0.5203 


1.226 

1.0215 


0.6379 

0.5315 


0.9719 
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